Giant Colloidal Diffusivity on Corrugated Optical Vortices 
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A single colloidal sphere circulating around a periodically modulated optical vortex trap can enter 
a dynamical state in which it intermittently alternates between freely running around the ring-like 
optical vortex and becoming trapped in local potential energy minima. Velocity fluctuations in this 
randomly switching state still are characterized by a linear Einstein-like diffusion law, but with an 
effective diffusion coefficient that is enhanced by more than two orders of magnitude. 

PACS numbers: 05.40.-a, 05.60.-k, 82.70.Dd, 87.80.Cc, 42.40.Jv 



Brownian particles moving on tilted washboard poten- 
tials exhibit two well-characterized limiting behaviors . 
When the potential energy wells are deeper than the ther- 
mal energy scale, diffusing particles become trapped in 
local minima. Their long-time self-diffusion coefficient 
vanishes in this limit. At the other extreme, tilting the 
washboard steeply enough to eliminate potential energy 
barriers allows particles run freely downhill. Because 
diffusion is decoupled from translation at low Reynolds 
numbers, a freely running particle exhibits displacement 
fluctuations characterized by its equilibrium self-diffusion 
coefficient, Dq. Between the trapped and running limits, 
particles intermittently switch between the two states, 
drifting downhill at a mean speed set by the rate at which 
particles are thermally activated over barriers. The tra- 
jectory, x(t), of such an intermittently trapped Brownian 
particle nevertheless is predicted 0, El 0] to satisfy the 
usual Einstein relation, 



lim (x 2 (t)) - (x(t)) 2 = 2Dt, 



(1) 



where D is an effective diffusion coefficient. Equally sur- 
prising is the prediction 0, 0, El D can be enhanced by 
orders of magnitude over Dq at the crossover from locked 
to running states. 

In this Letter, we provide experimental confirmation of 
substrate-mediated giant diffusivity by tracking the mo- 
tions of a single colloidal sphere traveling on tilted wash- 
board potentials created with corrugated optical vortex 
traps la, |6j, |jfl . The optically driven particle undergoes 
normal diffusion even in the intermittent regime, with an 
effective diffusion coefficient that increases more than a 
hundred-fold at the point of maximum intermittency. 

Our samples consist of colloidal polystyrene spheres 
2a = 1.48 jum in diameter (Bangs Laboratories, lot num- 
ber 6064) dispersed in water and confined within a glass 
sample volume formed by bonding a #1 coverslip to a 
microscope slide. This assembly is mounted on the stage 
of a Zeiss S100TV Axiovert inverted optical microscope 
for observation. Images are captured by an NEC TI- 
324II video camera and recorded on a Panasonic DMR- 
E100H digital video recorder for processing and analysis. 
The polystyrene spheres sediment into a layer roughly 
200 nm above the coverslip Individual particles are 



o 


Id) n 

\°) a = 




a = 0.1 


(b) ^ 


We 











3 
2 
1 



3 

2 5" 
1 5 
g 

3 
2 
1 





71 



[radians] 



FIG. 1: (a) Optical vortex with I = 40. (b) and (c) Cor- 
rugated optical vortices, a = 0.1 and 1.0, respectively, (d) 
Circumferential intensity profiles, 1(6), measured from (a), 
(b) and (c) at the radius of maximum intensity, Re. 



clearly resolved with a 100 x NA 1.4 SPlan-Apo oil im- 
mersion objective lens that also is used to project holo- 
graphic optical traps 0, H, 0, El into the sample. 
We tracked particles' motions with 10 nm spatial resolu- 
tion at 1/30 s intervals using standard methods of digi- 
tal video microscopy |l3|. From measurements on freely 
diffusing spheres, we estimate 0, Q| a wall-corrected 
self-diffusion coefficient of Dq = 0.19 ± 0.02 ^m 2 /s. 

Tilted washboard potentials were created from super- 
positions of ring-like optical traps known as optical vor- 
tices 0, 0, ■ Each optical vortex in this superpo- 
sition is formed from a helical mode of light whose 
fields, 

Mr) = u(r)e M , (2) 
are characterized by a radially symmetric amplitude pro- 
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file, u(r), and a phase f(r) = £9 proportional to the angle 
9 about the optical axis. The wavefronts of such a beam 
take the form of an £-io\d helix whose pitch determines 
the radius, Re, of the projected ring of light A 
typical optical vortex with £ = 40 and Rg = 4.2 ± 0.1 fim 
appears in Fig. ^a). The focused ring of light acts like 
an optical gradient force trap, drawing nearby dielectric 
particles to its circumference. The helical pitch also en- 
dows each photon in the beam with lh orbital angular 
momentum |l8j that can be transferred to an illuminated 
object. The resulting torque causes a trapped object to 
circulate around the ring |20| . 

Superposing optical vortices with opposite helicities, £ 
and —i, creates corrugated optical vortices such as the 
examples in Figs.^b) and (c) whose circumferential pro- 
files are sinusoidally modulated with 2£ intensity max- 
ima @ . An even superposition creates a so-called optical 
cogwheel [ilj consisting of bright spots arranged in a cir- 
cle of radius Re- This superposition carries no net orbital 
angular momentum and thus exerts no torque. A more 
general superposition, 



V>(r) 



V>i(r) + y/aijj-e(r) 
V / T+a 



(3) 



still exerts a torque, but also has a sinusoidal corrugation 
whose depth is set by < a < 1. The two limits, a = 
and a = 1, correspond to an ideal optical vortex and an 
optical cogwheel, respectively. We vary a by calculating 
phase-only holograms encoding the desired superposition 
using the direct search algorithm ^1 an d projecting the 
results with a Hamamatsu X7550 PAL-SLM spatial light 
modulator. Our experiments were performed with 1.5 W 
of laser light at a wavelength of A = 532 nm from a 
Coherent Verdi laser. 

The images of corrugated optical vortices in Fig. 
were captured by placing a mirror in the microscope's 
focal plane and capturing the reflected light with the ob- 
jective lens. Circumferential intensity profiles, 1(9) = 
\ip(Ri,9)\ measured from these images are plotted in 
Fig. ^d) and reveal errors in a smaller than 5 percent. 
Additional intensity variations of roughly 10 percent arise 
from imperfections in the optical train and so are in- 
dependent of a. These variations have an even smaller 
effect on the potential energy landscape experienced by 
the particle because the particle's finite extent tends to 
smooth them over |22j, |23( . 

A photon in the superposed beam has probability 
1/(1 + a) to have orbital angular momentum +£h and 
probability a/(l + a) to have orbital angular momentum 
—£%. The corrugated optical vortex therefore carries a 
local orbital angular momentum flux (£\/c) [(1 — a)/(l + 
a)] 1(0), where c is the speed of light. A fraction of this 
orbital angular momentum is transferred to a trapped 
object, and drives it around the ring. Circumferential in- 
tensity gradients, dgl(9), modulate this torque, and also 




FIG. 2: Instantaneous angular speed for a single 1.48 /im 
diameter polystyrene particle circulating around corrugated 
optical vortices at £ = 40, P = 1.5 W and a = 0, 0.36, 0.5 
and 1.0. 



induce optical gradient forces. The overall circumferen- 
tial force therefore has the form 

F(9) = F [A(a) + B(a) cos(2£6») + 77(61)] , (4) 



with A(a) 



(1 - 
+ a). 



a) /(l + a) and B(a) = 
The force scale, F is propor- 



2ya(A*+e)/(l 
tional to the laser beam's power and dep ends on the par- 
ticle's shape, size and composition [2j||23j. Given this, 
the washboard's tilt depends only on the excess, A(a) 
of right-helical photons in the corrugated optical vortex. 
The sinusoidal term, by contrast, includes a material- 
and geometry-dependent constant, £, accounting for the 
relative strength of the optical gradient force. We have 
omitted an irrelevant phase from the sinusoidal term's 
argument, and will ignore landscape's roughness, rj(9), 
in what follows. 

Figure |21 shows brief samples from typical single- 
particle trajectories in the trapped, running and intermit- 
tent regimes. Here, we have plotted the angular displace- 
ment, A9(t) = 9(t + St)-9(t) over the period St = 1/30 s 
of one video frame. The driving term vanishes in the cog- 
wheel limit, a = 1, and the particle remains trapped in 
a single local minimum of the potential, where it under- 
goes thermally driven fluctuations about its equilibrium 
position. In the freely running limit, a — 0, it circulates 
around the ring nearly three times a second. Periodic fea- 
tures in the running state's displacements result from the 
particle passing repeatedly over the disordered landscape, 
r}{6). From these, we estimate (^(9)^ ~ 0.01. At in- 
termediate values of a, the particle makes thermally acti- 
vated transitions between trapped and running states so 
that its trajectory is characterized by intermittent bursts 
of motion resembling random telegraph switching noise. 
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FIG. 3: (a) Dependence of circulation frequency on depth 
of corrugations. The solid curve is a fit to Eq. © yielding 
£ = 0.19 and Fq — 1.5 pN. The dashed curve is an estimate 
for uj(i(a) from Eq. @. (b) Diffusion coefficients estimated 
from the inset plots of mean-squared displacements (triangles) 
show a dramatic enhancement around a — 0.4. Squares de- 
note equivalent results obtained from Eq. JSJ . Long-term self- 
diffusion coefficients are unmeasurably small in the trapped 
state for a > 0.6. The solid curve is a guide to the eye. The 
dashed curve is a comparison to Eq. Q using parameters 
from (a). 



Averaging A9(t) over a period long compared with the 
trajectory's correlation times yields an estimate for the 
mean circulation rate u)(a), plotted in Fig. Eta). The 
net orbital angular momentum flux driving this circula- 
tion decreases as a increases, and we estimate the free 
circulation rate to be wrj(a) = A(a)u>(0) from Eq. |@J. 
This is drawn as a dashed curve in Fig. 02a) . The solid 
curve in Fig.|3ta) is a comparison to Stratonovich's exact 
expression for the mean drift velocity 0, , 



u)(a) 



-2tt 0R e F o A(a) 



where 



(5) 



H± P ±PV{9) 

Do 



and where V'{6) = F [B(a) cos(2^0) + T)(6)]. 

Slow instrumental drifts in u)(a) amounting to a few 
percent over several minutes can be estimated and 
subtracted off to reveal the linear growth of mean- 
squared positional fluctuations shown in the inset to 
Fig. Efb). The associated effective diffusion coefficients, 
D(a) are plotted as triangles in Fig. HJb). These values 
agree well with those obtained 0, Ej from fluctuations in 
the time, Tj(a), required to complete the j-th circuit, 



D(a) = 2ir 2 Rj 



(6) 



which are plotted as squares. 

The measured effective diffusion coefficient agrees with 
the equilibrium value, D(a) = Dq when the particle is 
in the free-running state a < 0.2. Larger values of a 
correspond to deeper corrugations that tend to trap the 
particle for longer periods. Longer periods of localiza- 
tion might be expected to reduce the particle's effective 
diffusion coefficient. Indeed D(a) — when the particle 
is trapped altogether for a > 0.6. Instead, intermittent 
trapping dramatically increases the effective diffusion co- 
efficient, with D exceeding 100 Dq at a = 0.4. 

This extraordinary substrate-mediated enhancement 
of the effective diffusivity is accounted for by the exact 
formulation, analogous to Eq. JSJ, due to Reimann et al. 

S0, 



D = D, 



11(6)1.(6) § 



(7) 



e T/3V(8T<t>) e -/3ReFoA(a) 



Analysis of Eq. (7J reveals that peak diffusivity should oc- 
cur when the driving force is just barely balanced by the 
periodic modulation. Neglecting r/(0), this occurs when 
a satisfies A(a) = B(a). Estimating the peak position to 
be a = 0.41 ± 0.01, we obtain £ = 0.19 ± 0.01. The sharp 
peak in D(a) therefore can be used to probe how ab- 
sorption and scattering transfer orbital angular momen- 
tum to objects trapped in optical vortices. The form for 
D(a) predicted by Eq. (J7J) is plotted as the dashed curve 
in Fig. Otb). Whereas Eq. (J5J agrees well with to(a), 
Eq. (JJJ describes a much broader and weaker peak in 
D(a) than is observed experimentally. It appears, there- 
fore, that D(a) is exceptionally sensitive to details of the 
tilted washboard potential. 

Histograms of angular displacements in Fig. 0] provide 
insights into this sensitivity. A freely running particle's 
displacements fall into a nearly Gaussian distribution, 
whose width and peak position both increase linearly 
in time. In the intermittent state, by contrast, parti- 
cles spend much of their time localized in traps, so that 
the short-time displacement probability is highly non- 
Gaussian, as shown in Fig. EJa). Because the poten- 
tial energy landscape is periodic and the effective par- 
ticle density is fixed, the displacement probability must 
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FIG. 4: Displacement histograms in the intermittent regime 
at q = 0.36. (a) Short time, t = 33 ms. (b) Long-time, 
t — (T) — 1.3 s. 



evolve into a Gaussian distribution through the central 
limit theorem. This self-averaging can be effective over 
as little as a single mean first-passage time (T) , as shown 
in Fig. Efb). This accounts for the essentially normal 
diffusion evident in the inset to Fig. |5ft)1 . 

The peaked structure of the probability distribution in 
Fig. Il^a) also suggests a qualitative explanation for the 
overall enhancement of diffusivity. A particle undergo- 
ing intermittent transport has a probability p(t) to be 
trapped for time t, and a probability I — p(t) to travel 
with an angular speed ujq set by the washboard's overall 
tilt. The mean drift speed in this highly simplified two- 
state model is oj = (I — p) ujq. The running state is char- 
acterized by thermal fluctuations in the mean-squared 
angular speed of magnitude Do/(R^t). This, however, 
can be dominated by fluctuations due to thermally ac- 
tivated transitions between the stationary and running 
states. Taking t = it/(£ujo) to be the time required to 
travel between potential wells in the running state, 



D 



Rft 
p(l-p) irR 2 e 



(<- 2 > 



ujq + (l-p) D . 



(8) 
(9) 



The effective diffusion coefficient therefore can be made 
arbitrarily large by increasing u>q, with an upper limit set 
by the onset of inertial effects. 

Giant diffusivity can degrade the performance of sort- 
ing methods such as gel electrophoresis and optical frac- 
tionation that exploit differential transport through a 
structured medium. In particular, the relative distribu- 
tion Ax/L of a sample that has traveled a mean dis- 
tance L through the landscape at speed v is y/2D/(Lv), 
which can diverge with the effective diffusion coefficient, 
D. This effect may be responsible for anomalous band 
broadening in electrochromatography [25| . Figure EU[b) 
demonstrates, however, that undesirable dispersal due to 
giant diffusivity can be overcome by more rapid driving, 
and that a small increase in driving force can have a dis- 
proportionately large effect on sorting resolution. On the 
other hand, substrate-mediated giant diffusivity should 



be useful for thoroughly mixing and dispersing materials 
in microfluidic environments, and might also provide a 
strategy for enhanced mixing in granular materials. 
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